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Abstract 

The dynamics of TeVeS in a homogeneous and isotropic universe is shown to be equivalent to the dynamics 
of an interacting two-component system, consisting of a scalar field and a "fluid" , related to the matter part, 
with explicitly given coupling term. Scaling solutions (solutions with a constant ratio of the energy densities 
of both components) in the "Einstein frame" are found which are exponentially expanding or contracting 
with no remaining freedom for Bekenstein's F function. In the "physical frame" these solutions are of the 
power-law type. An equivalent General Relativity (GR) picture of the dynamics suggests that it is the scalar 
field which plays the role of dark matter, while the "matter" has to mimic (phantom-type) dark energy. 
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INTRODUCTION 



The late-time cosmological acceleration, observationally suggested in has been backed 
up over the years by growing direct and indirect evidence. Direct evidence is provided by the 

a n 

luminosity-distance data of supernovae of type la (SNIa) [2j (see, however, [3j]), indirect evidence 
comes from the anisotropy spectrum of the cosmic microwave background radiation (CMBR) [4], 



from large-scale-structure data [5J, from the integrated Sachs- Wolfe effect [6], from baryonic acous- 



tic oscillations [7|] and from gravitational lensing 



. An interpretation of these results within Ein- 



stein's GR requires the existence of an unknown dark sector that contributes roughly 95% to the 
total cosmic energy density and manifests itself only through its gravitational action. Usually, this 
sector is divided into two components, dark matter (DM) and dark energy (DE), although there 
exist unified models of the cosmic substratum in which one single component shares propertied of 

n 

both DM and DE (see [9j and references therein). While the concept of DE is rather new, DM has 
been a topic both in cosmology and astrophysics for a long time. It was introduced by Zwicky as 
early as 1933 [lo| . In a cosmological context DM is required to explain cosmic structure formation. 
It has to provide the potential wells for the baryonic perturbations after radiation decoupling. By 
itself, the fractional baryonic energy density perturbations would not have had sufficient time to 



enter the non-linear regime. Despite of many attempts and the existence of severa" 
physical nature of DM remains a mystery until today (for recent reviews see 11 



candidates the 



13]) . In many 



occasions its gravitational action in GR-based cosmology is simulated by a pressureless fluid com- 
ponent. The failure of detecting DM over the decades has provoked more radical approaches which 
question the universal validity of GR. The hope here is, that within a (still to be found) modified 
theory of gravity DM may be superfluous. One remarkable step in this direction was Milgrom's 
MOND (Modified Newtonian Dynamics) [l4j. The intention of MOND is to have visible matter 
as the only source of gravity and to avoid the necessity of a dark matter component. The main 
feature of this theory is that it modifies Newtonian gravity for accelerations below a threshold of 
the order of 10~ 10 m/s 2 . MOND has been quite successful in describing galaxy rotation curves. 
For a recent review see [lj|. On the other hand, MOND is a non-relativistic theory and cannot 
account for a cosmological dynamics. A relativistic version of MOND, called TeVeS (Tensor Vector 

n 

Scalar), was developed by Bekenstein [161 ] . 

As are scalar-tensor theories, TeVeS is characterized by two different metrics. But while in 
scalar-tensor theories these metrics are related by a conformal transformation, the corresponding 
transformation in TeVeS is more complicated. In particular, it includes an additional unit time- 



2 



like vector field. Moreover, there appears a free function in TeVeS, that is not determined by the 
theory. Einstein's GR can be recovered from TeVeS in a suitable limit. Cosmological applications 
of TeVeS were studied in 



2014291] . In several of these studies the question was raised whether TeVeS 
can also account for an accelerated expansion or, in other words, whether it can also avoid the 
explicit introduction of a dark energy component [3, [3]. Reviews on TeVeS can be found, e.g., 
in [it], [ijj • A critical account of MOND-like gravity was given in [lj]] . 

Because of its many degrees of freedom TeVeS is technically complicated. Postulating for the 
physical metric a structure of the Robertson- Walker type with spatially flat sections simplifies the 
formalism but the resulting dynamics is still not straightforwardly transparent. It is the purpose of 
this paper to obtain simple analytic homogeneous and isotropic background solutions which allow 
for a discussion of at least some of the aspects of TeVeS. Our strategy here is similar to that already 
applied previously for scalar-tensor theories [js(3]. We shall make use of the circumstance that in 
the frame which is the analogue of the Einstein frame, the basic dynamics can be mapped onto 
the dynamics of an interacting two-component system. Then we look for scaling solutions of this 
system, i.e., for solutions with a constant ratio of the energy densities of both components. For 
previous studies of scaling solutions of the cosmological dynamics see, e.g., Subsequently, 
these solutions are transformed into the "physical" frame, resulting in power-law solutions for 
the physical scale factor. Already known solutions are rediscovered as special cases in a different 
context. Of course, power-law solution can only be considered as approximations for a specific 
epoch. In particular, they cannot describe any scenario, i.e., the transition between periods of 
different equations of state of the cosmic medium. Each of these phases has to be characterized by a 
different power-law solution, always under the condition that a power-law solutions is a reasonable 
approximation at all. Based on the explicitly known "physical" frame dynamics of TeVeS we 
establish an effective GR equivalent and compare the emerging picture with the predictions of the 
ACDM model. The perhaps surprising results of our calculations are: (i) an early matter phase can 
only be recovered in the presence of an additional phantom-type dark energy, although the latter 
is dynamically negligible, (ii) To be consistent with the observations, a separate phantom dark- 
energy component is required for the present Universe. In our solution TeVeS does not geometrize 
dark energy. But in accordance with the general intention of MOND type theories, no separate 
dark matter is necessary. TeVeS can geometrize cosmological dark matter, but only in the presence 
of an additional phantom dark energy component. While it remains open whether or not solutions 
of this type are (at least approximately) realized in our Universe, we hope that they may at least 
shed some light on the structure of TeVeS. 
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The paper is organized as follows. Section HH recalls the basic structure of TeVeS. In section ITTT1 
the general theory is applied to homogeneous and isotropic cosmology. In section IIVI we map the 
corresponding dynamics on the dynamics of an interacting two-component system with a coupling 
term that is explicitly prescribed by the structure of TeVeS. In section [V] we find scaling solutions 
for this system, i.e. solutions with a constant ratio of the energy densities of both components. An 
effective GR description of the cosmic medium is given in IVIl Section [VIII summarizes the results 
of the paper. 



II. GENERAL THEORY 



The dynamical quantities in Bekenstein's TeVeS are Einstein's metric tensor g^ u , a timelike 
unity four-vector U a and a scalar field (p. The physical metric tensor is constructed out of 
these quantities by 



24> , 



=207 



(1) 



The action of the system is 



S — Sn + S s + S v + S m , 



(2) 



where 



1 



9 16ttG 

with Newton's gravitational constant G, 



(3) 



- / d^x^g 



a 2 h^4>,^,u + ^F(kGa 2 ) 



(4) 



where a is a non-dynamical scalar field, A; is a dimensionless constant, I is a constant length 
scale and h^ v = g^ v — U IJ 'U V (the use of hf* v instead of g^ v avoids superluminal propagation of 
perturbations) , 



S v 



—ttG I 
32 J 



d x*J— ~g~ 



2A 



K 



(5) 



where K is another dimensionless constant and A is a Lagrangean multiplier. The square brackets 
in U\ a i £, denote antisymmetrization. Finally, S m is the matter part of the action, 



S m = j d A x^/^L r 



(6) 



4 



with the matter Lagrangean L m = L m (guv, •••)• "Matter" is understood here in a general sense. 
In particular, it can also be a dark energy component. The basic dynamic equations of TeVeS 
are obtained by varying the action ([2]) with respect to U^, a, (j), g^ u and A. In general form they 



are given, e.g., m 



19] . Here we are interested in the spatially homogeneous and isotropic case for 



which the basic set of equations considerably simplifies. 

III. BASICS OF HOMOGENEOUS AND ISOTROPIC COSMOLOGY 

The physical metric is assumed to be of the spatially flat FLRW type, 

ds 2 = -dP + a 2 [dr 2 + r 2 (d9 2 + sin 2 9d^ 2 )] . (7) 

The conversion relations between the quantities with tilde (physical frame) and without tilde 
(Einstein frame) are 

dt = exp y>] dt , a = exp [—(j)} a . (8) 

With the definitions 

da da 

H = — , H = — , (9) 
a a 

the relation between H, the Hubble rate in the Einstein frame and H, the Hubble rate in the 
physical frame, is 



H = exp[-cf>]\H--^) . (10) 

The matter part is assumed to be describable by the energy-momentum tensor of a perfect fluid, 

f a/3 = pu a U/3 + p (g a/ 3 + UaUp) . (11) 

At this stage it is left open, which type of matter T a p is supposed to describe. As already mentioned, 
it may also have a negative equation-of-state (EoS) parameter and account for some form of DE. 
The matter four- velocity and the vector Un are related by = e^U^ with = 5^ [if]]. The 
matter energy density is 

Pm = f 00 = p exp[20] , (12) 



the matter pressure 



p m = p exp[-20] . (13) 



Varying the action ([2]) with respect to g^ u results, when specified to the spatially homogeneous 
and isotropic case ([7]), in a Friedmann type equation for the Einstein frame scale factor a, 



1 fda\ 2 8vrG 



with an effective energy density 



and in 



p eff = p exp[-2<^] + a 2 (^j 2 + ^ F{kGo 2 ) (15) 



2d 2 a 1 f da\ 2 „ . , 



with an effective pressure 



Ccr 4 

p eff = p exp[-20] + a* ( J j - F(fcG^ 2 ) . (17) 



The physical matter energy density p obeys the equation 



\ (it a dt J 
which is the rewritten conservation equation 

| = -3#(p + p) (19) 

of the physical frame. The scalar field equation, obtained by the variation of §2$ with respect to 
0, specifies to 

Finally, varying ([2]) with respect to a yields 

G d , A) 



cr 



= h£ • (21) 



4/ 2 At 2 v ; \dt 

Equations JHJ), (fl~6j) . (|20|) and ([2T]) are the basic equations of TeVeS in a homogeneous and isotropic, 
spatially flat universe. 

IV. A TWO-COMPONENT DESCRIPTION 

The structures of the expressions (fT5|) and (fT?]) suggest a split (cf. {20! ]) 

p e // = 94, + Pf (22) 
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and 



Peff =P<t>+Pf , (23) 
respectively, with the field part (subscript (j>) 



and a fluid part (subscript /) 



as well as 



Pf = p exp[-20] = p m exp[-4</>] , (25) 



(26) 



x (it / 4Z 2 
and 

p/ = p exp[-20] = p m . (27) 

The fluid EoS parameter is 

P -l = P eXp[ " 2 ^ = g ^ t5 ■ (28) 
P/ /? exp[-20j) p 

We shall restrict ourselves here to constant values of the fluid EoS parameter w. The field EoS 
parameter is 

a = ; W , (29) 



(7 



We emphasize (cf. [2l|), that does not coincide with the matter energy density which is p m = Too 
(cf. Eq. ()12p ). Consequently, the decomposition (j22[) is not a split into field and matter parts. The 
theory does not admit a clear separation between the both. Also the fluid parts (I25D and (|28p do 
depend on 0. However, according to (EHJ) , the ratio g coincides with the physical EoS parameter 
w = S. As already mentioned, negative values of u> are admitted here. 

The coupling between the components can be made transparent by considering the correspond- 
ing energy balances. With the help of ([2~T]) the dynamics of the previous section is equivalent 
to 



and 

dp t 3 da . dd> , . , . 

The individual balance equations (I30p and (I31D add up to the total energy conservation 

dp e ff 3 da , , 



It is the system of equations (J30J) and (|3ip which we are going to study in the following. This 
system consists of two components which are interacting with each other. The interaction is given 
explicitly by the right-hand sides of equations (j30|) and (|3ip . It vanishes in the particular case 
-h. It is expedient to emphasize that equation (12ip (sometimes neglected in the literature 



w 



2 If ] ) was used in order to arrive at the system ([30]) and (f3Tj) . We shall solve the system of equations 
(|30p and (|3ip under the condition w = constant. Subsequently, we shall also assume = constant 
as well. Our aim then will be to look for scaling solutions, i.e. solutions, for which the energy 
density ratio y- is constant. 

V. A SCALING SOLUTION 
A. General properties 

For w = constant, equation (f3Tj) can be solved immediately, 

Pf = Pfo exp[(l + 3w)4>] I — J , (33) 



a 

where the subscript denotes some reference value. Via (|25p and (]12p corresponding expressions 
for p™ and p, respectively, are straightforwardly obtained. The latter reproduces the solution (24) 



m 



2l|. Defining a ratio R by 



P<t> 

the balance (|30]l becomes 



R = ?t , (34) 



pj, dt a dt dt 
For the simplest case that both and R are constants, equation (f35j) has the solution 

-(l+3iB)0fl /^0\ 3 ( 1+ ^) 



(36) 
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From now on we shall restrict our discussion to this case, i.e., to constant values of and 
R. While constant EoS parameters can certainly not account for a continuous transition between 
different epochs of the cosmological evolution such as the transition from decelerated to accelerated 
expansion, they are useful as simple exact solutions which are valid piecewise. Scaling solutions 
of the cosmological dynamics, i.e., solutions with a constant ratio of the energy densities of the 
dynamically relevant components, have been studied, e.g., in 



31 



38(| . We mention also that the 



ratio of the energy density of the scalar field to that of the matter, corresponding to R^ 1 in our 
notation, turned out to be approximately constant for Bekenstein's choice of the free function F 
(cf. [17]). This may be seen as an additional motivation to investigate such kind of solutions. 
It will turn out, however, that the requirement of a constant R does not reproduce Bekenstein's 
choice of F but leads to a different expression. Studying scaling solutions of the mentioned type 
also implies that we neglect the baryon component of the Universe since its contribution to the 
total energy budget is known to be small and the baryons as well as the also omitted photons are 
not expected to participate in the scaling of the dynamically dominating components. 
Combination of (j33j) and (f36l) yields 

Pfo (i+mmi+ji) (2o\ 3 ^ 



P<po ^ a / 

For this expression to be constant we may consider two cases. 

(i) The ratio R in ([3"T|) is constant for <fi = const and w = w^. Since <fi = const implies = — 1 
(cf. (f29l) ). we have also w = —1. Then p^, pj and p e ff are constant and a cx e ±Ht with H = const. 
Now we take into account the scalar field equation (|20p which can be written as 

An exact solution for = const requires pf = — > R = (unless w = —5), which means it is 
a vacuum solution. This solution reproduces the relations (22) and (23) in [211 ]. However, there 
exists another solution for pj = 0, namely a 2 = 0. As we shall see, this solution corresponds to 
Wcfi = —1. The point here is, that in the limit R — > it follows that a 1 — > 0, but, on the other 
hand, a 4 F — > const as will be shown below. Under these conditions (f29|) leads to = —1. 

(ii) The second case for which R in (137ft is constant is (we assume ao = 1 from now on) 

e (i+3«))0(i+R) = a 3(«i-tB^) ^ e (i+3w)<p = aS-r+rf _ ( 39 ) 

This relation fixes the dynamics of <fi. It implies 

# _ „ w-wt 

~dt ~ JTTMSjJTTT) ' (40) 



i.e., a proportionality between -& and H. Again we exclude here the case w = Using (|39p in 
the expressions ([36]) and (f33j) for p^ and pf, respectively, we find the scaling solution 

l+w +H(l+tg) 

P0 cc pf cc p e ff cc a i+« . (41) 

Via Friedmann's equation (114p it follows that 

ffaia^ i+« . (42) 

Notice that both signs are possible here. There are both expanding (H > 0) and contracting 
(H < 0) solutions. Integrating (|32"|). we obtain 

a<xt \ °ff) , w eff = ^ + • (43) 

w e // is the EoS parameter in the Einstein frame. Notice that it has the same structure as a 
corresponding term in scalar-tensor theories jsoj]. From (j39|) we have (to 7^ —4) 



— ef f "7 "\ 

a 3 ^W" => e ^ cx t ( 1 +-e//)d+3») _ ( 44 ) 



Further constraints can again be obtained from the scalar field equation (|20p . which is now written 

as 

Changing on the left-hand side the variables from i to a, we obtain 

1H^H (1+W) ^- (46) 

With H from (j4"2j) and p/ from (|4T|) . integration of yields 



3 



ff-r = — ( 1 + 3w))^ T - : -a - ' r+H (47) 



where pjo an d Hq are the values for pj and H respectively, for a = 1. Taking into account also 
(f42l and 3Hg = 8TiGp e ffo, relation (147]) can be transformed into 

1 - W0 + (1 - w) y 1 



Comparing the previously found relation (|40p between -S and H with ([48 p . allows us to obtain 



SttG 2 = g (l + ^)(l + 3u;) 2 

3 9 (i2 - u^) [1 - u> + 22 (1 - w)] " 1 J 
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Equation (|49p implies that the parameter a 2 is constant. Now, if a 2 is constant, the function F is 
constant as well. Via (|2ip . a constant F means that = const. But with F = constant and = 
constant, the energy density (|24p is constant, hence (j25j) and (j22j) are constant. This property 
leads to a condition on the solution ([IT]) namely, 

p e // cx a 1 + R = const = —1 — R (1 + w) . (50) 

Inserting (|50p into (|49p the latter becomes 

8vrG 2 i2 (1 + 3*) 2 

cr = 7 — " w — r • (51) 

3 18 (1 + to) (1 + R) K ' 

The quantity a 2 is related to the fraction and the EoS of the matter. It diverges for w = — 1. In 
the vacuum limit R = we recover the previously mentioned solution a 2 = (see the discussion 
following eq. ([38)) ). 

Now we consider the case w = — i for which the equations (|30p and (|3ip decouple. The solutions 
for pf and then are pf oc a -2 and oc a~ V + ™«v ; respectively. In such ft. CftSC ft constant energy 
density ratio requires = — | as well. But with to^ = tt; = — | the last relation in (|50p cannot 
be satisfied for a non-negative ratio R. Consequently, there does not exist a scaling solution of the 
type discussed here for w = This case will be excluded from now on. 

Use of the last relation of (j50j) in the expression (j4~3]) for w e ff, consistently yields w e ff = — 1 
together with H = constant (cf.Eq. (pEj) ). Then a oc exp[i7t], replacing (03]), and the relations 
dSI) reduce to 

= a 3 ^W oc e 3// ^* . (52) 

If we use (|50p in ()40p we obtain 

^ = 3±±^H. (53) 
dt l + 3w; v 1 



21| |. which was introduced there as an ansatz 



The first relation (|52p reproduces equation (25) in 
to simplify the dynamics without a further reasoning. Here it appears as a property of our scaling 
solution. For w > — | as well as for < — 1 the quantities and H have the same sign, for 
— 1 < w < — | the sign is different. Relation (|53p is valid for any constant value of R. 

From (|5ip we can draw the conclusion that a 2 is only positive for w < — 1. This is exactly 
the condition under which the authors of 2l| obtained exponential expansion. A negative a 2 , 
necessarily obtained for w > — 1, makes the kinetic term in (|24p negative, thus describing phantom 
energy. This is consistent with ([50|) . where < — 1 implies u) > — 1. In particular, this property 

11 



holds also also for a non-relativistic matter EoS w = 0. As we shall see shortly, despite of the 
negative kinetic term, the energy densities of the components and the total energy density are 
always positive. 

Phantom type solutions (within GR) are known to suffer from stability problems 



39 



40]. On 



the other hand, there are suggestions from higher dimensional theories that a phantom character 
"on the brane" may be an effective phenomenon 



41 



42J. To the best of our our knowledge 



phantom solutions for TeVeS and their properties have not been discussed so far. In any case, 
since a phantom component is not excluded observationally, this issue deserves attention. 

Another interesting special case is = 0, equivalent to w < —1. Under this condition both 
components change their roles compared with the case w = 0. Now, the scalar field mimics matter 
(although not separately conserved), while the fluid component is of the phantom type. Since 
TeVeS is a relativistic MOND-type theory and the intention of MOND was to make dark matter 
superfluous, it is now the scalar field itself which replaces DM. The "fluid" on the other hand, has 
to be identified with a DE component. We emphasize, however, that all these considerations refer 
to the dynamics in the Einstein frame. 



B. The equation of state 



For the total EoS we find, using (15011 . 

Peff _P$+Pf 



1 



> W <P + 



R 



-w 



1 



(54) 



Peff P<t> + Pf 1 + l + R 

The total EoS parameter is always —1, consistent with p e ff = const. Assuming the matter to 
behave non-relativistically, i.e., w = 0, corresponding to = — 1 — R (cf. (|29p). we find from (|5"Tj) 
and (1531) that 



8vrG 



1 



R 



-a 



dct> 

and — 
18 1 + R dt 



The corresponding relations for the special case 
8vrG 



-a 



i(i+M! 



and 



= 3H (w = 0). (55) 

0, equivalent to w = — 1 — 1/R, are 

H K = 0) • (56) 



1 



3 2 l + R " dt l + 2i?/3 

With R = one realizes the existence of a vacuum solution = const of Eq. (|30p also for = 0. 

Friedmann's equation (fT4"|) can be written 3H 2 = 8itG [p<p + Pf] with p^ and pj given by ([2 
and ([25]) . respectively. Then, with (f34"|) . it is equivalent to 



3H 2 = 8vrG (1 + R) 



dt 



(57) 
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Upon combining (|5ip and fj53[) we arrive at 



1 dt ) 8?rG 21 + iJ v ' 



for the first term in the bracket in (|57p . Introducing (|58p in (|57p leads to 



. r G 4 OEr2 l + f (l + t£» 

^ J F = 3g 1 + • (59) 



Combination of (1581) and (1591) results in 



Vrft/ l + #fl + *)4/ 2 V 7 



1 + f (1+tD)4J 
Consequently, the energy density ([21]) becomes 



1 ^ in 2 2 



i + #(i + ^^-i^TT+^UJ • (61) 

The expression ([6T]) . hence the total energy density ([22]) and the energy density ([25]) . are positive 
for any positive F and w > —1 as long as a 2 < 0. Consequently, although a 2 < 0, equivalent to 
a negative kinetic term, the energy density of each of the components as well as the total energy 
density are always positive. For w$ = 0, i.e. w < — 1, we have a 2 > and 

Pt/> = 2a 2 ( w<t> = 0) . (62) 



The pressure ([26]) becomes 



l + R(l + w)G 4 _ 1 + 12(1 + «;) 2 /#\ 



2 



l + #(l + *)4^ fTF = 2 *(! + «>) * UJ ' (63) 

together with (|6ip consistently reproducing the EoS parameter (150p . Thus, the Einstein frame 
dynamics of the system is completely solved for the given configuration. In the following subsec- 
tion we shall use the transformation relations ([8]) and (I10p to obtain the corresponding physical 
dynamics. 



C. Dynamics in the physical frame 

The transformations from Einstein's to the physical frame are mediated by exponentials of <f>. 
With ([52]) one has 4> = const for w = —1. For the physical scale factor, relations ([8]) and ([52]) 
provide us with 

a = e'^a => a = a !+ 3 * . (64) 
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For any * < — |, the scale factor a is proportional to a positive power of a. It is only for * = — 1 
that a = a. For non-relativistic matter we find 

1 1 

a = -p =>- a = -T7T (* = 0) . (65) 

In this case, a growing (decaying) a generates a decaying (growing) a. In particular, to have an 
increasing scale factor a, the corresponding a is the scale factor of a contracting solution. For 
Wfj, = 0, equivalent to w < — 1, we have 

j Wtsny a = 5( 1+3 /( 2fi )) (w = 0) . (66) 



1 

a = a 1 " 



A growing solution always remains a growing solution. Generally, we have p = p/exp [2(f)]. With 
pj = const (cf. ((HI) with ([50])), it follows that 

p oc <r 3 ( 1+ ^ . (67) 

This dependence is consistent with the balance (|19p . Notice that this is the dynamics of the matter 
component under the condition of a total EoS = —1, where the matter energy density p m 
via (I25p and (|22l) is part of the total energy density p e ff- It is not just the dynamics on a given 
background. 

In a next step we consider the relation (fTUj) between # and H. With (f53l) and ([52]) it follows 
that 

1 + 3* 1 + 3* K ' 

Recall that w ^ —1/3. For any w > —1/3 the rates H and H have different signs. Then, a 
contracting H corresponds to an expanding H and vice versa. The EoS w = — 1 is the only case 
for which H = H. For non-relativistic matter the connection between both Hubble rates is 

H = T = -2a 3/2 H (w = 0) . (69) 

For H = const < 0, the physical Hubble rate increases with a 3 / 2 . This demonstrates again the 
phantom character of the solution. On the other hand, for = 0, equivalent to a phantom-type 
EoS of the matter, we obtain 

H = L__ a - 1 /( 1+2fi / 3 ) H K = 0) • (70) 

l + 3/(2i2) \ <P ) \ ) 

In this case, H and H have always the same sign. An explicit expression of a in terms of t is 
obtained with the help of (|68p with (|64p . Then, 

H = 2 —~ a l^)H => —„ = 2 —~ a \^)R (71) 

1 + 3* dt 1 + 3* v ; 
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with a constant H. The integration yields 



2 

3(l+ii) 



(72) 



We mention that a solution a oc t 3 ( 1 +™) has been obtained before in [21[. According to (|7ip . an 
expanding solution for the physical Hubble rate, H > 0, corresponds to -ff < for w > — g. But 
with H < and u> > — | the denominator in (|72p vanishes for t > ti. The solution (|72p describes 
an expanding universe that approaches a singularity after a finite time (big rip). The time span 
is determined by the vanishing of the denominator in (|72p . For the case of non-relativistic matter 
f2l) reduces to 



1 + 35? H (t-U) 



(w = 0) 



(73) 



The big rip singularity is approached after a time span ij, r 



1 + 35? H (t br ~ U) = 



1 



tbr ti ~\~ 

m\H\ 



(74) 



Notice that this big rip solution is different from Caldwell's solution in GR [43J. There, it is a 
consequence of w < —1. In our case it follows for w > — ^ which implies a contracting solution 
H < in the Einstein frame. No singularity occurs for w < — 1. For = 0, i.e., w < — 1 and 
i7 > 0, we find 

-i 2i?/3 



a = a,- 



1 + 



1 



H = °) • 



(75) 



1 + 2.R/3 

Accelerated expansion requires R > 3/2. 

The time coordinates t and t are related by the first transformation of ([8|). With (|52p we find 



3-H 



1+m 



t oc e n-3 



l + 3ru 

Ft oc lnts i+* 



(76) 



as well as 



e v oc - 



a 

a oc — 



Combining a oc exp [fl"t] with (|76p and (|77p , it follows consistently that 



a oc 



p(i+» e// ) 



(77) 



(78) 
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where we have introduced the physical EoS parameter 

w eff = -2-w. (79) 

For any w > — 1, including the case w = 0, the total equation of state is of the phantom type. For 
w = — 1 it follows w e fj = —1. It is only for w = —2 that we recover a total matter-type behavior 



eff = 0) which, according to ([50]) . implies w$ = — 1 + R. A radiation behavior requires w 



7 



For = 0, i.e., w < — 1, the parameter w e // becomes 

*e// = -l + ^ K = 0). (80) 

The medium as a whole can be seen as characterized by a physical energy density 

p eff oc a" 3 ( 1+ ^/) , (81) 

as solution of the conservation equation 

^ + 3H(l+w eff )p eff = 0. (82) 

With (|79p . the dynamics of the physical frame is explicitly known. The cosmic substratum as 
a whole behaves like a fluid with an effective EoS w e jf. The relation between the physical EoS 
parameter w e fj and the auxiliary quantity w e ff = — 1 is different from a corresponding relation 
in scalar-tensor theories [30]. In scalar tensor theories, both frames are related by a conformal 
transformation which is not the case in TeVeS. 

The deceleration parameter in the physical frame is 



IF V H dt ) \ L H dt J ^ W dt 



da \ ^ / 1 J_ dfji' " 

~dt) y H dt 

It is related to the corresponding parameter q = —-^jp^r in the Einstein frame by 

J_# / 1 dH i\ 1 d ( 1 d(j> 

H dt \ H dt ] \W dt 1 ) H dt \ H dt 



(83) 



1 _ 2 

1 H dt 



(84) 



We have accelerated expansion if the numerator of the right-hand side of eq. (|83ft is positive 



/ 2 dd>\ ( 1 dd>\ 1 dH 1 d 2 6 

For the dynamics discussed in this paper (w = const, = const, R = const) we have q = —1. 
With = ^| = in eq. ([83]), use of eq. §3§ provides us with 

? = -^ = ^(l + 30e//) • (86) 
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Consistently, we have q = — 1 for w = —1. For non-relativistic matter q = — | is valid. In this case 
q does not depend on R. For = 0, i.e. w < — 1, the deceleration parameter becomes 

9 = -(l-^) K = 0). (87) 

Accelerated expansion is obtained for R > | in accordance with (f75j) . The deceleration parameter 
vanishes for R = | which corresponds to w e // = ~~ § (cf- Eq. (|80p ). For i? = 1, equivalent to 
w = —2, we consistently recover q = ^. 



D. Implications for the F function 

Now let us combine Eqs. (I2ip and (153p . The result is 

o" y l + 3w 2 



(1 + 3*) 

Then, the relations ([88]) and (f59|) provide us with 



i^l _ 24, G i^SM; . . (89) 

a'F (l + 3to) 2 [1 + |(1+ _ 



W 



Using the expression (|5ip in (|89|) . the latter takes the form 



da 2 



(a A F) R l + w 



a 4 F 2a 2 1 + | (1 + w) 

with the solution 



(90) 



4.+3{l+w)R 

F oc (a 2 ) 2 +( 1+ *)« . (91) 
For the special case w = with <r 2 < we have 

3i?.+4 

F oc (-cr 2 ) «+ 2 («3 = 0) . (92) 

For Wcf) = 0, equivalent to w < — 1, the dependence of F is 

F oc ^ K = 0) . (93) 

This is consistent with = in (|29p under the condition ^ = const. 

While in general there does not exist a theory for Bekenstein's F function, we found that for 
the special situation considered in this paper, there is no remaining freedom for the choice of 
this function. The existence of a scaling solution requires a structure that is different from that 
suggested by Bekenstein. It remains open, whether the dependence (|91|) is a large scale limit of 
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a more general function. We mention that in the vacuum limit R = the relation (|92p implies 
a F = constant which proves the statement following eq. (j38[) . This is, together with a 2 = for 
R = 0, consistent with an EoS w$ = — 1 (cf. eqs. ([29]) and ([50]) ). It is also interesting to realize 
that the structure of (I92p is similar to that of an additional (to Bekenstein's ansatz for F) term 
introduced ad hoc in [20j]. For the special case R = the resulting inverse quadratic dependence 
on — a 2 in (j92p even exactly (up to the sign) coincides with the structure of that term (cf. eq. (21) 
in Q). 

VI. THE COSMIC MEDIUM 



The energy density (|8ip can be seen as the solution of a conservation equation (|82|) with an 
effective pressure 

Pes = w eS p eS , (94) 



p cS oc a 3 ^ . (95) 



where w e s is given by ([79]) . Then 



The overall dynamics is completely determined by the EoS parameter w of the matter. Does this 
solution have any relevance for our real Universe? As already mentioned, for any w > — 1 the 
effective energy density p e s increases with a, i.e., the total dynamics is of the phantom type. If, on 
the other hand, w itself is of the phantom type, i.e., w < —1, one has w e s > — 1 for the medium as 
a whole. The matter dynamics is described by (cf. Eq. (j!9p ) 

poca- 3{1+ ^ . (96) 

Equation (|82p with the solution (j8ip and (|79p represents a GR equivalent for the TeVeS dynamics. 
Given the separately conserved matter component with energy density p and the conserved total 
energy with a density p e ff > one may associate a so far unknown component with an energy density 
p y to the difference between p e g and p, 

PcS = P + P y ■ (97) 

Let us further assume that the geometry based y-component can be described by an effective (not 
necessarily constant) EoS parameter w y such that 

^l + (l + w y )p y = 0. (98) 
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In an effective GR description, the Universe is now made up of two separately conserved components 
with energy densities p and p y . Then 



Introducing the ratio 



we may write 



Since 



w e g = w-^— + Wy^r- . (99) 
PcS PcS 



f=-§~, (100) 



Py 



ui y = (1 + f) w e Q — fw . (101) 



p 

Pcff 



and 1 + f = ^ , (102) 



Pes-P l ~TT~ 

Pcff Pcff 



and, since the ratio J- is known, 

Peff 



we obtain 



For w y we find 



4- = n a-^ l+ ^ with n = -P- , (103) 

PcS PcSO 



= _^o« and x ~ = _ (1Q4) 

l-Ooa- 6 ^) 1 - n a- 6 ( 1+ ^ 



w y = w eS - If (1 + w) = -2 - w - If (1 + w) . (105) 
In terms of the scale-factor, 

2 + w ( i + n d-^ 1+ ^) 

(106) 



y 1 - ft a- 6(1+,5) 

As already mentioned, equations of state u> e // = § and u) e // = can be reproduced by tj = — | 
and a) = —2, respectively. For iu e // = | relation (|106p specifies to 

. 1 1 + 7n a 8 1 
o 1 — SZoO 05 
For a <C 1, i.e., at high redshifts, the EoS parameter w y approaches w y = |. The y component, 
an entity of geometric origin, behaves as radiation as well. With respect to the DM problem, the 
case w e ff = is of particular interest. For this case Eq. ()106j) takes the form 

= - _ | (™e// = 0) . (108) 
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For a <C 1 the EoS parameter w y approaches w y = 0. The y component behaves as non-relativistic 
matter. So, the apparently strange result is, that to have a non-relativistic matter period w e ff = 
we need a y component with w y = that only exists in the presence of a phantom type DE 
component with w = — 2. The total energy density is the sum of the energy densities of these 
components. With w = —2 the EoS of the scalar field component is = — 1 + R. It is only a 
combination of uu and R which is fixed. The apparently most natural choice = is realized 
for R = 1. But for the physical ratio f we obtain 

from (I104p . i.e., f — > for a <C 1. This is equivalent to 

Pj, > p => Peff-Py (5<1), (110) 

as one expects for a matter-dominated universe. The DE component p is dynamically negligible 
at high redshifts, but its presence is crucial for the non-relativistic matter component to exist. 
These results demonstrate, that physical equations of state such as w e ff = \ or w e ff = can 
be obtained within TeVeS on a purely geometric basis, provided only that there exists a separate 
phantom-type DE component. Finally we also mention that for w = —1 which implies w e tt = 1 
we obtain w y = — 1 as well. 



At the present epoch (I106P reduces to 



2 + tBfl + n 



™y0 = ~ • (HI) 

l — S'O 

Since we have shown, that there exists a solution with w y = at high redshifts, it seems tempting 
to consider the possibility of a present value w y o ~ as well. This corresponds to the original 
intention of TeVeS to avoid the introduction of a DM component. Under the condition w y o ~ we 
have from Eq. (jllip . 

2 , , 

w ps =- . (112) 

l + ^o 

If the y-component is to play the role of DM today, the "matter" component with an EoS w 
has again to be associated with DEJn an equivalent GR picture. As already mentioned, this is 
in accordance with the findings in 



21| where the authors obtained an inflationary solution for 
w < — 1. Now we know that the actual background universe is quite well described by the ACDM 
model. Within the latter, the total EoS parameter is 

^ACDM = _ (n3) 

£// 1 + PM0/PA 
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where ~ |. This corresponds to a value Jlo ~ f- Using this ratio in (I112p . we find an 
effective EoS for the DE equivalent of w ~ —I = —1.14. Such a value is compatible with current 

nn 

observational data [44J, [45[ . Under the condition (|112|) it follows from ([50]) that 

+ + - . (114) 

It is only a combination of to^ and R which is fixed, not each of these quantities by itself. However, 
the value for w is consistent with = which requires R = 7. With this choice it is obvious that 
the scalar field represents the DM. 

These considerations show that different epochs of the cosmological evolution are reproduced 
for different values of w. An early radiation period requires w = — |, a subsequent matter period 
has w = — 2 and for the present epoch we found tt; ~ — y. This indicates a tendency to less negative 
values of w with the cosmic time. Guided by the ACDM model, one may speculate that in the 
future w approaches w = —1, which means w e ff = — 1 as well. As in the ACDM model, the 
final state of the universe would be of the de Sitter type. As already mentioned, the parameter 
w y , according to (|1U5|) . also tends to w y = — 1. In this limit, all the components, including 
the y— component, behave as DE. The final state of the Universe would be characterized by two 
different DE contributions, one from the matter sector, the other one of geometric origin. However, 
since power-law solutions do not allow to construct a scenario, this remains speculative. Perhaps 
a dynamical system analysis could clarify whether there exist fixed points with corresponding 
properties. On the other hand, any value of w different from w = — 1 in the long-time limit would 
invalidate the split (I97p . Namely, for any w < — 1 the ratio jyj in (I103p would grow to values 
larger than unity, which implies a negative energy density p y together with a negative value of 
f in (|104p . Consequently, a limit w = — 1 is simply necessary for our power-law solution to be 
applicable for the future universe with a > 1. 



VII. DISCUSSION 



We have solved a simple case of the TeVeS dynamics in a homogeneous, isotropic and spatially 
flat universe. The total effective energy density was split into two components in a way that allowed 
us to look for scaling solutions, i.e., solutions with a constant ratio of the energy densities of these 
components. One of the components is a scalar field with a potential part, related to Bekenstein's F 
function, the other one represents a "fluid" contribution. TeVeS dynamics in the "Einstein frame" 
then corresponds to an interacting two-component dynamics with explicitly given coupling term. 
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The Einstein frame Hubble parameter is necessarily constant and Bekenstein's F function is fixed. 
For the physical Hubble parameter we find a power-law solution, corresponding to a total effective 
energy density p e ff with an EoS parameter w e ff = —2—w, where w is the "matter" EoS parameter. 
"Matter" is understood here in a general sense which includes w < as well. We establish an 
effective GR picture of the dynamics according to which the cosmic substratum consists of two 
separately conserved components: the first one is the mentioned "matter" with energy density 
p and EoS parameter w, the second one is characterized by an energy density p y = p e ff — p of 
geometrical origin and a (generally not constant) EoS parameter w y which is calculated from the 
explicitly known dynamics of p e ff and p. A non-relativistic matter period of the cosmic expansion 
can only be obtained for w = —2. The "matter" has to represent (phantom type) dark energy. 
For small values of the scale factor we find a non-relativistic matter- type phase w e ff = with 
w y — > 0, i.e., the geometry based component acts as matter. Moreover, the dark energy fluid is 
dynamically irrelevant in this limit. However, its mere existence is absolutely necessary to obtain 
a realistic dynamics. This demonstrates that an equation of state for non-relativistic matter can 
be obtained within TeVeS on a purely geometrical basis. For the actual Universe the situation is 
less clear. But, of course, it has to be approximated by a different power-law solution. Assuming 
tentatively the ACDM motivated present value ~ ^° - ~ | together with w y o = 0, this requires a 
present dark energy EoS parameter of w ~ —1.14, which seems not yet ruled out by the data. 
In a sense, this picture, although perhaps surprising, is consistent with the original intention of 
MOND type theories, to avoid the introduction of a dark matter component. While our solution 
implies a geometrization of dark matter, the existence of a separate dark energy component is 
not only necessary to obtain the observed amount of accelerated expansion but also to recover a 
non-relativistic matter phase of the cosmic medium. Future investigations along these lines will 
have to include baryon and radiation components as well as perturbations about the homogeneous 
and isotropic background. 
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